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Abstract: Adding fundamental matter of mass mq to = 4 Yang Mills theory, we 
study quarkonium, and "generalized quarkonium" containing light adjoint particles. At 
large 't Hooft coupling the states of spin < 1 are anomalously light (Kruczenski et al., 
hep-th/0304032). We examine their form factors, and show these hadrons are unlike any 
known in QCD. By a traditional yardstick they appear infinite in size (as with strings in flat 
space) but we show that this is a failure of the yardstick. All of the hadrons are actually of 
finite size ~ g'^N /mq , regardless of their radial excitation level and of how many valence 
adjoint particles they contain. Certain form factors for spin-1 quarkonia vanish in the large- 
g^N limit; thus these hadrons resemble neither the observed J quarkonium states nor p 
mesons. 
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1. Introduction 

The discovery that gauge theory at large 't Hooft coupling can be described using string 
theory has given rise to many interesting developments in both fields. The original 
idea has been extended to purely four-dimensional confining gauge theories Q ^ These 
examples provide the first toy models of QCD whose kinematics is that of the real world. 
In particular, these theories share with QCD the property of four-dimensional approximate 
scale invariance in the ultraviolet, as well as confinement in the infrared. Consequently, 
the differences between these models and QCD are largely due to dynamical, rather than 
kinematical, issues. Indeed, these theories smoothly become QCD-like when the 't Hooft 
coupling is small. 

However, most of the theories studied up to now differ from QCD in a significant way: 
their matter comes in representations with of order N'^ fields, and they are all neutral under 
some portion of the center of the gauge group. These differences from QCD are much more 
important, both kinematically and dynamically, than the fact that many of these theories 
are supersymmetric. The reasons for this are clear. Quarks in the N representation of 
SU{N) color introduce new features into the 1/A^ expansion; quark pair-production makes 
confining flux tubes unstable; baryons now appear in the spectrum; and there are new flavor 
symmetries which may or may not be explicitly or spontaneously broken. Thus the theory 
without quarks in the N representation, supersymmetric or not, differs signflcantly from the 
theory which contains them. It is thus important to the development of these toy models of 
QCD to introduce matter in the fundamental representation. In the limit where the number 
of flavors Nf is much less than N, this was considered in ||^, ^, 0; a simpler method was 



invented in and studied further in 10|. Still more recently, there has been additional 



work in other related contexts [11|. 



One of the obvious objects to study in a theory with quarks is quarkonium. At small 
't Hooft coupling this is just a hydrogenic atom, but at large 't Hooft coupling the system 
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is highly relativistic. The quarkonium spectrum at large 't Hooft coupling, in the model 
studied in [^] , is remarkable from the field-theoretic point of view. Its details could not 
have been guessed from any known theoretical argument or from any aspect of the observed 
QCD spectrum, and it is profoundly tied up with the representation of the four-dimensional 
gauge theory as a higher-dimensional string theory. Still more puzzles emerge in states built 
from a quark, an antiquark, and one or more particles of the much lighter adjoint matter. 
These states would naively be expected to be qualitatively different in size and structure 
from pure quarkonium. Instead, it has been found that these states are quite similar 
to pure quarkonium at large 't Hooft coupling. 

In order to gain better insight into the structure and couplings of these bound states, 
we have computed some of their form factors and transition matrix elements with respect 
to various conserved currents. Fourier transforming the form factors, we find that all the 
hadrons are more or less the same size — larger than one would expect for quarkonium, 
and smaller than one would expect when adjoint matter is bound to the quark-antiquark 
system. Our results do not solve any mysteries, but they do raise interesting questions 
and suggest other calculations to do in future. We are also led to a conjecture about the 
substructure of the hadrons. 



We review the results of Q and |10| in the following section. Sections 3, 4 and 5 contain 
our methodology, a summary of our results, and the detailed computations. The physical 
implications of our results are discussed in Section 6. 



2. Preliminaries 

The introduction of matter in the fundamental representation into theories with gravita- 
tional dual descriptions has been considered by a number of authors ^, ^. However, 
theoretical prejudices about the appropriate systems, and technical difficulties with those 
that were investigated, delayed progress for some time. Recently, Karch and Katz cut 
the Gordion knot, pointing out that many interesting questions could be addressed in the 
simplest possible brane construction with fundamental matter: a small number A'^^ of D7 
branes in the vicinity of a large number N of D3 branes. 



2.1 The theory in question 

The field theory corresponding to this arrangement of branes consists of = 4 SU{N) 
Yang-Mills coupled in an AA = 2 supersymmetric fashion to Nf hypermultiplets in the 
fundamental representation. We will write the = 4 vector multiplet as an AA = 1 
vector multiplet Wa and three chiral multiplets ^1,^2,^3 in the adjoint represntation. 
The hypermultiplets are given as AA = 1 chiral multiplets Q^,Qr (r = 1, . . . ,Nf) in the 
fundamental and antifundamental representation respectively, and we will call their scalars 
"squarks" and "antisquarks" and their fermions "quarks" and "antiquarks." Written in 
M = 1 language, the theory consists of kinetic terms for all the fields (along with M = 1 
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superpartners of the kinetic terms) and a superpotential of the form 

Nf 

W = V2 tr([$i, «>2]$3) + J2 Q"'^3Qr + mrQ'Qr 

r=l 

where is the mass of hypermultiplet r and the trace is over color indices. 

The theory has an 50(4) ~ SU{2) x SU (2) symmetry, consisting of an SU{2)^ symme- 
try rotating $i and $2 and an SU{2)r J\f = 2 R-symmetry. The charges of the fields under 
these symmetries are shown in the table, where we write = + iX^ , <I>2 = + iX"^ ■, 
and <I>3 = + iX^, the adjoint fermions as Ai, A2, A3, A4, and the quarks as "0; V'- 





5C/(2)$ X SU{2)n 


U{1)r 


U{Nf) 









1 




(0,0) 


±2 


1 


Ai, A2 


(io) 


-1 


1 


A3, A4 


(0,^) 


+1 


1 




(0,^) 





Nf 




(0,0) 


Tl 





The masses are the eigenvalues of a matrix which transforms in the adjoint of U{Nf). 
If all masses rrir are zero, the U{\)b. M = 1 R-symmetry is preserved by the superpotential. 
Nonzero masses break this U{1) R-symmetry, but leave invariant the 50(4). The U{Nf) 
is also generally broken by the masses, but if all the masses are equal, as we will assume 
throughout, U{Nf) is preserved. 

The string dual description of this field theory ||8|, for large 't Hooft coupling,^ is simply 
given as Nf D7 branes placed as probes inside AdS^ x S^. Each D7 brane introduces a 
hypermultiplet of mass m, where m is proportional to the minimum value of the AdS^ 
radial coordinate r to which the D7 brane descends. At any larger and fixed value of r, 
the D7 brane fills out an subspace of the S^. We will give more details on this brane 
construction below. 

The theory's reduced supersymmetry is of little concern, but it does seem at first to 
have a serious dynamical problem. The gauge coupling now has a positive beta function, 
making the ultraviolet definition of the theory problematic, and potentially destroying any 
hope of making sense of the theory using gauge-string duality. In particular, the dilaton 
is no longer constant. Where the dilaton (equivalently, the gauge coupling of the gauge 
theory) becomes of order 1, one ought to do an S-duality transformation; but this would 
turn the D7-branes into magnetic 7-branes, which are very difficult to handle. 

However, the gauge beta function is very small. In particular, the beta function for 
qym is of order gyM^f^ ^° beta function for the 't Hooft coupling A is of order \^Nf/N. 

^We use the following notation: qym is the A/" = 4 gauge coupling, a = Qym/^''^^ ^ = the 't Hooft 

coupling, Qa is the type TIB closed string coupling, and the couplings are related under gauge/string duality 
by gs = a. 
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Another way to say this is that at large N and small A'^^ this is a naturally "quenched" 
theory: all quantum effects due to the fundamental matter are parametrically suppressed 
by Nf/N. The squarks and quarks are simply probes of the dynamics of the = 4 sector. 
From the supergravity point of view, this suppression is essentially the statement that the 
backreaction, through the variation of the dilaton, metric and other fields, caused by the 
D7-branes is negligible out to exponentially large AdS radius. Consequently, as Karch and 
Katz argued, we can study to leading order in Nf/N all of the issues that can normally be 
investigated in quenched QCD, including the spectrum of hadrons. 

2.2 The (s)quark-anti(s) quark hadrons 

Karch, Katz and Weiner |^] considered the formation of heavy-light mesons in this theory, 
working at large 't Hooft coupling A ^ 1. They found a number of surprising results, in- 
cluding confinement of quarks without formation of flux tubes (what one might call "Gribov 
confinement." |l2|) 

Another natural issue for study is quarkonium, namely bound states of a quark and 
antiquark of equal mass. At weak coupling (i.e., small A) this is simply the problem of the 
hydrogen atom, but at strong coupling we do not a priori know what to expect. 

A seemingly different problem is that of a quark and antiquark bound to one or more 
of the $i particles. In the simplest D3-D7 system, these particles are massless. An easier 
system to think about is obtained by giving the small bare masses, so that (as in M = 
confinement of flux occurs at a distance very long compared to the quarkonium mass 
and length scales. Then we can imagine stable bound states of a quark, antiquark and 
some number of the much lighter $j particles. This problem has not been much studied at 
small 't Hooft coupling, but one would expect these bound states, in which the quark and 
antiquark combined must be in the adjoint representation, would be very different in size 
and shape from the simple quarkonium bound states. 

Surprisingly, at large 't Hooft coupling both problems can be addressed simultaneously, 
and on the same footing. This was done by Kruczenski, Mateos, Myers and Winters pO[| . 
They derived the spectrum and higher-dimensional wave functions of mesons consisting 
of one (s)quark and one anti(s)quark of equal mass, along with some number of massless 
adjoint particles. The spectrum for states of large spin {s » \/A) was found to be the 
same as that of hydrogen; this is to be expected, since in this limit the (s)quarks move 
nonrelativistically even though the coupling is strong. For states of lower spin it was shown 
that the spectrum has a Regge-like relation betweeen mass and spin: m? oc s. And states of 
spin 1, ^ and were found to be extraordinarily deeply bound — so deep that the mass of 
the hadron is a tiny fraction of the mass of its (s)quark constituents. It is these anomalously 
light states of low spin — states that are described in the string theory variables by modes 
of massless higher-dimensional fields — that are the focus of this article. We now review 
their properties in detail. 

We generally follow the setup in . The near-horizon geometry of D3-branes filling 
the 0123 directions of the ten-dimensional space, and placed at the origin of the 456789 
coordinates, is given by AdS^ x S^. We write the metric on the Poincare' patch variously 
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9 

ds"" = ■^{dxi'f + ^dr^ + R^d^l = L^rj^^dxV + ^ ^{dx^f 

c=4 

where c = 4,5,6,7,8,9, = 0,1,2,3, = Ec(^'')^ and i?^ = a'^/^T^gJ^. Note r = is 
the horizon of and r = oo is its boundary. The D7-branes fill the 01234567 directions; 
each is placed at some position x^ + ix^ = rnqoi . We will make all the masses equal in this 
paper, so without loss of generality we can take the masses to be real (x^ = mqa' = L, 
x^ = 0.) The induced metric on the D7-branes is then 

ds"^ = -^■q^ydx^dx'' + Y_, 

c=A 

= ^{e^ + l)ry^,dx'^dx- + -^dg^ + R^^—dQl , (2.1) 

where = ^ — 1, and the involves the angular coordinates in the four-dimensional 
space spanned by 

We will use various coordinates for our calculations which we summarize here: 

V = {L/rf , w = l-v; g"^ = v'^ - 1 = . (2.2) 

Note the boundary of the space is located at r = oo, p = oo, v = and w = 1, whereas 
the point on the D7-brane which is closest to the AdS horizon at r = is at r = L, p = 0, 
V = 1 and w = 0. 

It is useful to define an overall mass scale 



rrih 



Barred quantities will be defined relative to rrih — for instance, we will often use dimen- 
sionless momenta q = ■\f\f^\lvnh and dimensionless masses fa = m/rrih. 

There are two real scalar fields and one gauge field on the D7 brane worldvolume; 
theses are the massless modes of the open strings whose ends are on the D7-branes.^ If 
there are Nf D7-branes, there are Nj such fields, transforming in the adjoint representation 
of U{Nf). The modes of the gauge field, when reduced to five dimensions, break up into 
spin-one and spin-zero modes in five dimensions. The spin-one modes are associated with 
the U(Nf) fiavor current, and operators built by adding $ fields to the current: 

{Q^^^Q - Q^^Q^)e0 = Q^^^di'Q + ^^^V'^Vq - Q^^^'^Q^ + • • • (2.3) 

where stands for any product of $i and $2 which is a symmetric and traceless represen- 
tation under 5*0 (4). These operators have conformal dimension Ai = £ + 3 and transform 
as (|, I) under 50(4) « SU{2)^ x SU{2)r, and as adjoint under U{Nf). 

^We use — h ++ signature. 

^Note these 7-7 strings are not dynamical in the gauge theory! In the Maldacena hmit of the D3-D7 
system, the 3-3 strings (the J\f = 4 sector) and the 3-7 strings (the hypermultiplets) are dynamical, whereas 
the closed strings and 7-7 strings act as background fields. In the dual string description, the 3-3 and 3-7 
strings are absent, and the closed and 7-7 strings are the dynamical degrees of freedom. 
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When these operators act on the vacuum they can create a spin-one bound state of a 
Q, a Q, and i <I> particles (plus any number of gluons, gluino pairs, and <&~anti-<^ pairs, of 
course.)^ In particular, for ^ = this is a spin-one quarkonium state, much the same as 
the J/^ or the T in QCD. The bound state spectrum, and the bulk wave functions of the 
individual bound states, were computed in where they were termed "type //". These 
hadrons are modes of the 8-dimensional gauge boson Am, of the form 

Ap = 0, ^« = 0, A^ = C^.ct>'\p)e'^-^y'{S^), k.C = (2.4) 



cPg = {Cli/R^)e\l + £>2)-i— ^F(-n, -l-n-tA + 2- -q") 
= (C/^/i?2)u>^/2(i _ w)-^''^F{-l - n-£,2 + n + e;i + 2;w) 

= (C'/„Vi22)^;(^+2)/2(l _ y)mpiW+l)(^2v - 1) , 

where n > 0, £ > 0, 



c£ ^,/2(2„.2.. 3) ('"-^^-^)(":'r^)^(":^r 



e + i j\ £ + 1 J \ £ + 1 y^^"' 

Here Pn"'^^ (x) denotes a Jacobi polynomial. The masses of these states are 

Mfj = Amlin + £ + l)(n + £ + 2) , 
where again mh = L/E? = mQ{X)^^/'^ . By coincidence the £ = states also appear in the 



model of |13|. 



Another mode which is easy to identify is the mode /— in [|10[, whose conformal di- 
mension is Aq = £ + 1 {£ > 1), and which transforms as (^y'-, ^-y') under 5*0(4) ~ 
5C/(2)<j, X SU{2)ji. For ^ = 1 it has dimension 2 and is a triplet under SU{2)r. This 
uniquely picks it out as a multiplet containing the M = 1 chiral operator 

m'-'Q)e,g=o = Q'^'''Q + --- (2-5) 

along with other modes related by 5*0(4) (mainly given by replacing Q with , $i with 
<I>2 or etc.) The mode /— is a mode of the 8-dimensional gauge field with its component 
on the 5^ nonzero: 

A^ = 0, A, = 0, A^ =ct>'-{p)e^^--y'^-{S^) (2.6) 



C = {CijL)e'+\l + g')-^-^-'F{-n,l-n-£-£ + 2--e') 

= {CljL)w'-^ (1 - u;)^(^~3)/2F(1 -n-£,2 + n + £-£ + 2-w) 

= {CljL)v^'+^y\l - vf+'y^Pjf+'''-'\2v - 1) 

^They cannot, in this theory, create a spin-zero state; this is not a general feature of the large 't Hooft 
coupling limit, however. 
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where n > 0, ^ > 1, and 



--^^<— <';:f)(";:t>(";::>^- 

These scalars have masses 

M]_ = Amlin + i){n + £ + 1) . 

Since we will not use the type /+ hadrons, we will henceforth discard the minus sign in I— 
when referring to these states. 

These are the two classes of (s)quark-containing operators and states that we will use 
for calculations in the rest of this paper. We will also consider matrix elements of the 
flavor current, which is the operator appearing in Eq. ( p. 3D with £ = 0. This will require us 
to know the nonnormalizable mode of the corresponding five-dimensional gauge boson at 
spacelike = (fiTif^, which is of the form 

4>Z^(w-a)= } ' Fi-a,l + a;2;w), 2.7 
sin(7raj 



where a = ^ (—1 + y/l — q^ 



For completeness we briefly comment about the other classes of hadrons appearing 
in |1C]. There is a set of complex operators + • • • which are obtained from the 

operators (|2.5|) by the action of two supersymmetry generators; these complex modes create 



the two "scalar modes" in [10[. The type /// and type /+ hadrons are harder to identify, 
as there are a number of candidate operators, of which only particular linear combinations 
are chiral.^ 

2.3 The = 4 sector 

We will also consider certain modes of the AA = 4 sector, which appear as states of the 
ten-dimensional supergravity on AdS^ x S^. We will only need the spin-one states created 
by acting on the vacuum with the 50(4) current J^^ ~ tr ^idf_i^j{T^Q^^-j)ij + • • • and the 

spin-two states given by acting on the vacuum with T^'^ ~ tr FpF^'^ + ^jd^d'^^ + • • • . (The 
traces are over color indices.) However, it will be convenient to consider also a spin zero 
operators which has a close relationship with the spin one and two operators. The spin 
zero operator corresponds to the ten-dimensional metric element which is a singlet under 
the transformation on the seven brane. Under the breaking SO (6) S'0(4) x 5*0(2), 
the traceless symmetric rank two tensor representation 20 branches to (1,1) © (9,1) © 
(1,2) © (4,2). The singlet corresponds to the dimension-two S'0(4)-singlet operator S ~ 



^Note there are no chiral operators of protected dimension containing "1>3 and Q together, because the 
form of the superpotential imphes c&aQ and Q^a are simply proportional to Q and Q, respectively, within 
the chiral ring. 
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2.3.1 The continuous spectrum of the conformal theory 

We will consider the operators whose spins range from zero to two. The non-normalizable 



mode for the operator with spin S is |14| 



where q = \f(^ jmh and > is spacelike. Note that 1? jF^v in the prefactor is the warp 
factor of the four dimensional part of the metric. It is needed in order to have the correct 
boundary behavior of the non-normalizable mode. 

The SO (4) current is associated under AdS/CFT duality with a five-dimensional gauge 
boson A, which descends from the off-diagonal elements of the ten-dimensional metric. 
In particular, let v"" be a Killing vector on the three-sphere which points purely parallel 
to the D7-brane world- volume (and thus leaves it invariant.) Then the metric elements 
h^a = A^Va,hra = ArVa (where ds"^ = ds\^g^^^g5 + h^adx^^dx"-) define the 50(4) gauge 
boson on AdS^. In radial gauge Aj- = 0, the non-normalizable modes of A at spacelike q^ 
are 

A^{q) = e^q^Kiiq^) (2.9) 



which are normalized following ( |2.8D so that A^ — > as u ^ 0. The corresponding 
normalizable modes, at specific values of timelike q^ = — \q^\ = —mfjh?, are of the form 

A^{q) oc e^fhy/vJi{rh\/v) . (2-10) 

One can build normalizable wave packets from the normalizable modes, although, just as 
with plane waves in flat space, there is no normalized version of the modes in ( 2.1C| ). 



For the spin-two current we need the modes of the traceless AdS^ part of the metric. 
We consider traceless fluctuations ds^ = ds^^ds^xs'^ ~^ h^ydx^dx'^ , = 0, in radial gauge 
h^^ = 0, hy^ = 0. The normalizable modes take the form 

h^y{q) oc e^yJ2{my/v) (2.11) 

while the nonnormalizable modes are 

KM) = ^f^vj^ ■ ^K2iqV^) . (2.12) 

As — > 0, the non-normalizable mode approaches hf^i, — > e^u^:^} which is the four- 
dimensional warp factor, times a polarization tensor. 

Since the = 4 theory is conformal in the infrared, its spectrum of normalizable modes 
is continuous: any timelike q^ is allowed. We should not think of the normalizable modes 
as hadrons, but simply as a spectral decomposition of a conformal field. 

2.3.2 The discrete spectrum of hadrons of a confining model 

However, it will often be useful to consider instead a confining theory, in which the conformal 
invariance of the AA = 4 sector is broken in the deep infrared. For instance, we can imagine 



-9- 



the = 4 sector is deformed into a theory similar to AA = 1* Q, where the N = \ chiral 
superfields $i,$2)^3 are given a small mass. This causes confinement to set in at a low 
scale A. In all known models of this type, the AdS radial coordinate is effectively cut off at 
r = r^m = Ai?^. Our results will be insensitive to the details of this cutoff, as will become 
clear, so a crude model will suffice to capture the essence of the physics. In any such model, 
the flavorless N = A sector of the theory will have a discrete spectrum of spin-s (s < 2) 
hadrons, with masses of order nvrA (0 < n E Z), and with mode functions given at large r 
by the normalizable modes of the corresponding bulk flelds. 

To be definite, we will model confinement through the boundary condition that the ten- 
dimensional wave function of ecich lia-dron satisfies the Neumann condition at t — i^min-j 

that 

is, at Vmax = ^ ■ 111 this simple model, the wave functions for the hadrons created by 
the 50(4) current are precisely those of ( |2.1(]| ) for v < Vmax, with a quantization condition 
that only qf^ satisfying {^Vmax Ji {'rn^Vmax ) ) ' = 0, where again m = \/[^/m/i, are allowed. 
Thus we have a countably infinite series of modes, with mass = fh^Vmax = Co;n, > 1, 
where Co;n is the n*^ zero of Jq. The normalization is obtained from |15] 

/ dvd^ne^^^g^'''AMAM = l, (2.13) 
Jo 

where 

ds^ = e^^rj^ydx^dx"" + ds\. 



Thus we have a spectrum of states with normalized wave functions 



. I . V 1\JvlVmaxJ\{C,i^:n\JvlVmax) o-3^ V2v Ji{mn^/v) 

^^^'^-^='^ WMW:) = ^ MC0;n) ' ^'-''^ 



and masses 

mn = = CO;n ' -, ~ 7 • \^-^^) 

nih nih n>i V 4/ mh 

We will need the nonnormalizable modes satisfying the same boundary condition, which 

A) 



are (for spacelike = q^mf ^ 



A,iq) = e,qV^ !^K,iqV^) + /i(g-V^)^^| . (2.16) 
Similarly, the spin-two hadrons created by T^^'^ have wave functions given by 



J2{Cl;n\/v/Vmax) _ A\/2J2(m„^) 1 '7\ 



with masses 

m„ ^ A / 1\ A , , 

rhn = — = Cl;n ^ + 7 ^ • ^.18 

ruh rrih n»i \ 4 J nih 

where Ci-n are the zeroes of Ji. The normalization constant is obtained in a similar way as 
in Eq. (pTl3|) : 

/ dvd^ne''^^g''''g''^hMphNQ = l . (2.19) 
Jo 
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The nonnormalizable modes are now 

3. Methodology 

3.1 Definitions and Notation 

We will concern ourselves mainly with the form factors from the hadronic matrix elements 
of the U{Nf) flavor current J^, the 50(4) « SU{2)^ x SU{2)r current J^, and the energy- 
momentum tensor T'^'^. We will also consider matrix elements of the S'0(4)-singlet spin-zero 
dimension-two operator 5 oc $f + $| — 2$|, which is part of the conformal A/" = 4 sector. 

We will compute a large number of form factors, so it is important that our notation 
be clear. When the maximal information needs to be displayed, we will use the following 
notation for a form factor. Suppose we have a spin-zero (type I) or a spin-one (type //) 
initial state with £, rii quantum numbers. In our calculations the final state will always be 
of the same type and share the same i, though n2 / ni in general. The operator whose 
matrix element we are computing will be referred to by an index S = 0, 1,2 for the spin- 
0,1,2 operators S, J^, and T^^^, and by an index / for the flavor current J^. The max;imal 
notation will therefore be 

4!l (3.1) 

In many cases one or more indices will be clear from context and will be omitted. Finally, 
in computations involving spin-one initial and final states, there can be more than one form 
factor. We will label these with an obvious subscript; in this case, the indices S,i,ni,n2 
will never be needed, and will be left implicit. 

For scalar hadrons, the matrix elements of the spin-one SO{A) current can be written 

{n2,i,p'\J^mm,e,p) = (77'^- - q''q''/q'')ip + p\Fi%^{q^) . (3.2) 

A similar expression holds for the flavor current. 

Spin-one hadrons probed by a spin-one current have (in a Lorentz-invariant parity- 
conserving theory) three form factors, one each for electric monopole, magnetic dipole, and 
electric quadrupole couplings to the spin-one current. We label these with subscripts e, rn, q. 
In general the current matrix element (n2, C'l-^l* l^ij ^j^*) C); of spin-one mesons takes 
the form: 

i{ [(C • C)ip' + P)u - ip' ■ OC. - iP ■ CKuW - q^q''/q^)Fe{'f) + 

+ [{q-Oe-{q<)C]Fra{f)+ (3.3) 

+^[{P ■ C'W ■ C)iP+p).M'' - q^q''/q^)F,{if) } x 

x{2'Kf5\'£p,), 

i 

where e, and are the polarization vectors of the current, in-state, and out-state re- 
spectively, and Fe, F^, and Fq are electric, magnetic, and quadrupole form factors, with 
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all other indices suppressed. For reasons to be discussed later, in the supergravity limit 
Fe = Fjn and = 0, so our notation will usually be to drop the subscript except when 
necessary, writing Fni^2 the 50(4) electric form factor and replacing 1 with / for the 
flavor current. 

The formulas for the spin-two current are similar, although more complicated. Again 
we will find there is only one non-vanishing form factor both for spin-zero and spin-one 
hadrons, which we will label ini,n2 when necessary. 



3.2 Determining the shape of the hadrons 

Initially we will compute the form factors in momentum space, where they are Lorentz- 
invariant functions of (f'. To obtain information that is easier to interpret intuitively, we 
would like to reexpress the form factors in position space. There are ambiguities in how 
this is to be done, and problems which might arise at timelike where there are poles. 

A four-dimensional Fourier transform of the form factor has the feature that it is Lorentz 
invariant. However this function does not have a well-known physical interpretation. More- 
over one must consider large timelike where certain difficulties with the supergravity 
approximation will arise. 

A three-dimensional Fourier transform of the form factor at spacelike is useful for 
nonrelativistic systems. For a two-body nonrelativistic bound state, this quantity is the 
square of the wave function, and a similarly simple interpretation applies for many body 
system: for a spin-one current, it gives the three-dimensional distribution of the correspond- 
ing charge. But our bound states are highly relativistic (since their binding energy is so 
large), and (for small ^, n) their form factors are large even when y^lg^l is of order the hadron 
mass. We are therefore not confident that this interpretation extends straightforwardly to 
our case. 

By contrast, the two-dimensional Fourier transform of F{q'^) (for spacelike q^) does 
have an interpretation which is applicable for relativistic systems. For a hadron moving 
with extremely high momentum in the z direction, and choosing q>^ = (0, gl, 0) in the x — y 
plane, we define 

nxl) = e"^^-^^F{q^) 

This function is the hadron's two-dimensional "transverse charge distribution" p±{x±) in 
the x—y plane, times 2tt in our conventions. The usefulness of this interpretation stems from 



its connection with generalized parton distributions |16, 17, [l^], as shown with considerable 



care and rigor in |19]. This applies in our case even though the (s)quarks have large masses, 
as long as the hadron is ultrarelativistic compared to the quark mass scale. ^ 

In some cases we will find elegant closed- form expressions for F(x'^). In others we can 
still compute many of their general properties: the large- and small-g^ behavior of F{q'^), 



®Since the form factors are functions of only, one can convert from any one of these transforms to any 
other, at least when restricting to spacelike q^. Our choice of the two-dimensional transform is therefore 
somewhat arbitrary; however it gives unambiguous and interpretable information about the structure of the 
hadrons. 
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the large- and small-x^ behavior of its Fourier transform, and some characteristic measures 
of hadron shape and size. 

One classic measure of the size of a hadron is given by the second moment (r^) of the 
transverse charge distribution. For any given form-factor, we can compute this moment 
using 

{r')^A^^F{q') (3.4) 

(where the 4 replaces the often-used 6 because we are measuring a two-dimensional charge 
distribution.) However, this measure is not unique. Another measure is (r) itself, which 
cannot so easily be obtained from F{q^)] it is best extracted directly from the Fourier 
transform F{x\). 

3.3 Calculational techniques 

The required calculational techniques are well-established and straightforward; see for ex- 



ample [20|. Each hadron we will consider is a mode of a particular five-dimensional field, 
which itself is a mode of an eight-dimensional field on the D7 brane or of a ten-dimensional 
field in the bulk. The hadrons containing a (s)quark and anti(s)quark will be of the so-called 
type / or type // class described in the previous section, both of which descend from the 
gauge bosons on the D7-brane. To compute the matrix element of a current, we need to 
examine the five-dimensional field whose boundary value couples to that current. For Jj, 
this is the mode of the gauge boson on the D7-brane which has its index in spacetime and 
is constant on the . (Acting on the vacuum it creates the type // hadrons Eq. (^^) with 
£ = 0.) For J4 the gauge boson in question is the dimensional reduction of a ten-dimensional 
supergravity mode; as in [^] it is associated with a Killing vector on the <Z . For T^'^ 
we need the five-dimensional massless graviton. To compute the matrix element {out\0\in) 
requires knowlege of the trilinear interaction between the three five-dimensional modes cor- 
responding to the initial hadron, the final hadron, and the operator O. This interaction can 
be derived from the Born-Infeld action on the D7 brane. In the supergravity limit (large 
A), all of the interactions that we will require are obtained from the single term 

1 / dSx^/^^9^«e-^FM7vFpQ , (3.5) 

where D, the dilaton, will not play a role below and will be dropped from future equa- 
tions. Here M, A^, P, Q are curved 8-dimensional indices, the metric is induced from the 
ten-dimensional metric, and gs is the eight dimensional Yang-Mills coupling, y/gla' . 

The matrix element is then given by plugging into the appropriate trilinear vertex a non- 
normalizable mode of the field corresponding to the current, and the wave functions of 
the incoming and outgoing hadrons. Performing the integral over the eight dimensions, 
we obtain the answer required. In general, the integral over the and the integral over 
Minkowski space will be elementary. The important integral will be that over the AdS 
radial dimension, and will take the form (for spin-zero incoming and outgoing hadrons) 

/dv 
(l>non{q,v) (l)*outiv) (t>in{v) (3.6) 
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where we remind the reader that v = {L/r)"^. For spin-one hadrons the integrand is different 
only in small details. 



4. Summary of our main results 

In this section we summarize our results, and in the following section we present the detailed 
computations. 

4.1 Form Factors 

4.1.1 General results for all theories 

We begin with some observations which apparently follow from conformal invariance and 
large- AT alone. We have not derived them from general arguments, but it should be possible, 
and would be interesting, to do so. In particular, it is not yet known whether the following 
properties are true only in the supergravity limit and thus apply only at large 't Hooft 
coupling. 

The results below apply to any hadrons with the property that (1) their mass scale is 
large compared to any other scale which breaks conformal symmetry, (2) the mass scale 
which sets their masses breaks conformal symmetry only at order 1/N. Such hadrons will 
have wave functions which solve an equation in the background of a conformal theory. If 
the conformal sector has conserved currents in addition to the energy-momentum tensor, 
then we find the leading form factors of the energy-momentum tensor are related to those 
of the currents, which are in turn related to those of operators of spin zero and dimension 
two. In our case, for scalar hadrons, the form factors for S, J^,T^^^ satisfy 

F''\^) = ^^{if^F^'\e) (4.1) 

for S = 0,1,2. Taking a two-dimensional Fourier transform of these form factors gives 
functions of (where x = mhX±,) satisfying 

^^'^(x^) = |^[(x^)^^^^°^(x^)] • (4.2) 

These relations follow only from properties of the modes corresponding to the operators in 
the conformal sector, so they also apply when the incoming and outgoing hadrons have spin 
one. Note that 



x=oo 
x=0 



Since F^^^ is not badly divergent^ at x — ^ 0, this shows it must fall, at large x, as 1/x^ for 
ni = n2, and faster for ni ^ n2. One can see then that F^^^ falls as (x)^^(^+'^\ or faster; 
the bound is saturated for m = n2, and also sometimes holds for m ^ n2, S > 0. 



7^(o)<!(^^2-j _^ g 2{Ao 1) jg^j.gg which implies is logarithmically divergent for Ao = 2 and finite 
for larger Aq. 
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In contrast, the flavor current does not couple to the conformal sector, and is sensitive 
to the masses of the quarks. Its form factors fall off exponentially at large x. Conformal 
invariance at large requires oc F^^^, oc but supergravity imposes a stronger 

condition, namely that the two form factors are actually equal at large q (when they are 
both normalized to 1 at = 0.) 

In general, spin-one hadrons can have three form factors under spin-one currents: elec- 
tric monopole, magnetic dipole, and electric quadrupole. We find (similarly to that the 
anomalous magnetic dipole {Fm — Fg) and the quadrupole form factor (Fq) are zero in the 
large-A limit. The nonvanishing form factor Fg = F^ has the same large-x behavior as in 
the spin-zero case, though its large q^ behavior is This is true for both SO(4) and 

flavor currents, and follows from properties of the couplings between hadrons and currents 
in the supergravity limit. This will be discussed elsewhere |21]. 

4.1.2 General results for this theory 

In our particular theory, considerable simplification is obtained from the fact that the wave 
functions in the bulk are given by powers of v and (l—v) times Jacobi polynomials. The form 
factors are easily analyzed in position space using the two-dimensional Fourier transform, 
where they are given by integrals of the form (|3.6| ). 
In general, we find 

Fi'X.i^') oc ni(x2) log + n2(x2) , (4.3) 

where Hi and 112 are polynomials. The first, up to normalization, is simply given by the 
wave functions themselves, or derivatives thereof: 




(4.4) 



is of degree 2{£ + n), and, importantly, begins at order (x'^)^"^. The second polynomial 
n2(x^) is of degree 2£ — 1 and begins at order 1. It can be determined by noting that at 
large x all terms which grow faster than x^^^^^'^^ must cancel. We will determine 112 in 
some specific examples. 

At large x, the coefficient of the potentially-leading x^^^^^'^^ term is 

/ 4m4n2 ■ (4-5) 

Jo V 

The above integral, for S* = 0, is the normalization integral for the normalizable wave 
functions, and it thus vanishes for ni ^ n2- For other S it need only vanish for |ni — n2| > S, 
as follows from properties of the Jacobi polynomials. 

Some of these facts have natural momentum-space counterparts. The behavior at small 
X — in particular the absence of a logarithm multiplying x^-^ for j < £ — 1 — is associated 
with the requirement of conformal invariance that the momentum space form factor fall as 
Xjq^^ . Similarly, the (x^)"^'^"'"^) behavior of F^^'^ at large x follows from the fact that the 
expansion of F{q^) near zero is as a polynomial plus {q^)^ \ogq. 
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For spin-one hadrons, Eq. ( [4. 3D is unmodified, but Eq. (|4.4| ) becomes 



Hi 



2 



9^ 



9f 



,1-5 



(4.6) 



The Hi polynomial now starts at order (x^)^"*"^ to account for the corresponding change in 
the large (f' behavior, which is now \/q^^^^. 

It is useful to evaluate the 50(4) form factors at x = 0: 



^ ni,ni 



2(2ni + 2£ + 2±l) 



(4.7) 



where the upper (lower) sign applies to hadrons of spin one (zero) . 

The flavor form factors are less amenable to such a description due to their mathematical 
complexity. We do not have general results beyond those of the previous section, except 
that all such form factors can be written as a sum over a finite number of spin-one £ = 
hadron poles: 

1, J. qz _|_ 



n=0 



where ttt.^ = 4(n -|- l)(n -|- 2). The — (-I-) sign in the upper limit of the summation applies 
for spin-zero (spin-one) external hadron states. In position space this form factor can be 
written as a corresponding sum of KQ{mnx) Bessel functions. 

4.1.3 Ground state form factors for each £ 



Since P, 



1, the ground states have form factors proportional to 

^1 dv 



(4.9) 



This implies 







F^^'i^') = ni(x^) log i^jf^j + n2(x^ 
ni = (-i)^-Mc/of (xY-Hi + xY^^ 

and 112 is proportional to the polynomial V-\ defined in Eq. ( ^.25| ). From this we can 

00 



X 



(4.10) 



obtain Fqq^^ , 5 = 1,2, as described above. 



For spin-one, the electric form factor has a similar form, with 

ni = (-i)^-MQVF (x2)^+\i + xY+i 

and 112 is proportional to the polynomial defined in Eq. ( p. 25 ). 



(4.11) 



For the flavor case, c^qq in Eq. (4^) is given in Eq. (|5.q ) for spin zero hadrons and in 
Eq. (5.45) for spin one hadrons. The small-x behavior of F^^^ is similar to that of F^^\ 
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4.1.4 Diagonal form factors at large t 

Remarkably, although these formulas become complicated at large an alternative and 
surprisingly simple formula can be used instead. We find that for both spin zero and spin 
one hadrons, 

Fo?^^(x) ~ i ^ (4.12) 

for £ ^ 1. This misses the (x^)^^^logx term, but this is a negligible error at large I. In 
momentum space 

F^^'it) « ^^{q/V2fKs{q/^) . (4.13) 



For the flavor case, the large i limit of q q is given in Eqs. (5.7) and ( ^.46 ) for spin-zero 
and spin-one hadrons respectively. 

4.1.5 Diagonal form factors at large n 

At large n, a formula can be obtained which is valid for sufficiently large x. 

= Is ( 5 ) (" » ("^) 

This formula applies for both spin-zero and spin-one hadrons. While inaccurate for very 
small X, these formulas are normalizable for > 1, and we find empirically that this formula 
gives reliable answers for moments such as (x) and (x^). The apparent divergence at x = 
is not present for finite n (except for Type /, £ = 1) as is clear from Eq. ( [4.71) . 

For the flavor case, c^^^^^^ for both spin zero and spin one hadrons is given in Eq. (5.8). 



4.1.6 Results for some ofF-diagonal matrix elements 

The conformal invariance of the = 4 sector and the derivative relations between the Jacobi 
polynomials imply some additional relations for off-diagonal matrix elements between the 
ground state and an excited state at a given i. 



From this expression the form factors for 5 = 1,2 can be obtained from our general result 
(|4.2|). The same result holds for spin-one hadrons with replaced with C^^ . 



4.2 The transverse sizes of the hadrons 

Armed with this information we can compute certain moments, in particular (r) = m^^(x) 
and (r^) = m^^(x^), with respect to the flavor, 5*0(4) and energy-momentum distributions 
of these hadrons; we use subscripts /, 1, 2 for moments of the corresponding form factors 
pif) ^ F^^\ F^'^\ Where these moments are infinite we regulate them using finite A (see 
Sec. p. 3. 2 ); where they are finite we set A — > 0. 
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4.2.1 The ground states n = at general i 

For every we can compute properties of the ground state (n = 0). For the spin-zero 
states, we find 

(n) = 2(r2) 



TT V{1 + l)r(2£ + 2) _ 

5\ ' 



imh r(^)r(2£ + |) 



1 



8m2 ^ + 1 ' 
0.496... , , / V / X 0.697... 

where i?^ = 1 + 1/2 + 1/3 + ••• + !/£ is the £-th harmonic number. 
For the spin-one states, the formulas are similar: 

, , , ^ r(^ + f)r(2^ + 4) 



2mft r(^ + 2)r(2^ + 



8m| ' 

0.672... , ^ , , , , 0.697. 

~ {rf)i=o < {rf) < {rf)e=oo ~ 

ruh ruh 

4.2.2 The ground states at large I 

For n = 0, ^ — ^ oo, for both spin-zero and spin-one hadrons, 

in) = 2(r2) - 



{rl) 



1 



- 18 - 



{rf) 



0.697... 

rrih 

log 2 



f ml ■ 

4.2.3 The small-^ states at large n 

For both spin-zero and spin-one hadrons, the hmit n — > oo at smah and fixed i gives 

in) = 2{t2) — ; 



{rl) ■ ' 



8ml ^ 



0.61 



, 2(1 -log 2) 

V// ^ 2 • 

ml 

5. Computations 

In this section we derive the results outlined in the previous section. 
5.1 Spin zero hadrons 

We now proceed to study the spin-zero (type I) hadrons created by the operators Q^^~^Q 
(and other members of the same 50(4) multiplet.) 

5.1.1 Flavor current 

In order to study the (s)quark-anti(s)quark hadrons using the matrix elements of the flavor 
current, we need to consider a situation with A'^^ > 1. The hadrons in question transform 
in the adjoint of U{Nf), so for Nj = 1 they are neutral under the flavor current. We will 
consider A'^^ hypermultiplets of equal mass, which leaves the U{Nf) unbroken. The modes 
described in [|lO[, which only depend on the quadratic terms in the D7-brane gauge fields, 
are unchanged for A'^^ > 1, except for transforming under a nontrivial representation of 
flavor. However, the cubic interactions among the D7-brane gauge bosons give the main 
contribution to the flavor-current matrix elements. 

For the coupling of the flavor current to two spin-zero hadrons, the important term in 
the D7-brane Born-Infeld action is 
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where / is the structure constant of the group SU{Nf) and we have rescaled the gauge 
fields A^^ — > ggA^^ to obtain this form. In this form, the boundary (r = oo) value of the 
non-normalizable mode times the coupling gs is set to unity. 

Substituting the non-normalizable mode for A^, and normalizable modes for the mesons, 
we obtain the matrix element: 



The above integral can be done by partial integration, using the modes in Eq. ( p.6| ) and 
Eq. ( |2.7| ). The normalizable mode is a polynomial in w and J dwF{a -\- l,b + 1; c + l;w) = 
-^F{a,b; c;w). By integrating the non-normalizable mode and differentiating the product 
of normalizable modes repeatedly, this integral can be evaluated. The general form of the 
form factor obtained in this way is 

ni+n2+2e-l e 

where m,^ = 4(n -|- l){n + 2) is the squared mass of the n-th vector meson with £ = 0, and 
Cn,ni,n2 is a Constant independent of q. 

For the general case, we can get the large q'^ behavior in the following way: Remem- 
bering that we evaluate the integral by repeated partial integration, the ^-times integration 
of 4)l^^{w, a) gives a term with its leading behavior 4^£!/(g^)^ for large if' and the function 
we need to differentiate repeatedly has the structure 

plus corrections suppressed by additional factors of (1 — w). Therefore, the first i — 1 terms 
obtained by partial integrations vanish and the leading large if behavior is given by 



) .2^^2{ni+i){ni+£ + l){n2+£){n2+e+l){ii^Y (g2)Ao-i' ^^""'^ 

where Aq = £ -|- 1 is the conformal dimension of the spin-0 mode. This behavior agrees 
with the "quark counting rules" in QCD [22, |2^, because the latter follow from conformal 
invariance and are not in fact limited to weak coupling and the valence-parton model of 
hadrons. 

It is possible to compute the Fourier transformation of the form factor. Using the 2-d 
Fourier transformation 

( 2^ 2 ) = ^oM, 
\q + m'^ J 
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we get 



ni+n2+2£-l 



mi, 



Yj cl,ni,n2^0(mnX) 



n=0 



2 

711 C 

h^h,ni,n2 



TT 



2mnX 



n=0 



fc=0 



(5.4) 



ni+n2+2f-l oo i^ m.nx 'i^fc 



where n is the smallest value of n with nonvanishing c^^.^^ ^^- If X]c^,ni,n2 ~ 0, the log- 
arithmic divergence for x — > vanishes; as can be seen from Eq. ( |5.2D , this is also the 
condition that the coefficient of at large vanishes. We know this must occur for 
i > 1. In general, from the fact that the form factor falls as l/q^^, we can deduce that 
S m n2("^")^"' = for j = 0, 1, • • • ,i — 2, and thus the leading non-analytic term in the 
small-x expansion is — [x^*^^"^) + • • • ] logx. 

In principle, we can get all c^^„^ by partial integrations, but it's difficult to get a 
closed form for „^ ™ general. In the case of ni = n2 = 0, we can get a closed form^ 
which is given by 



'^n,0,0 



Co!nC[o^B{i + l,n + e + 2)3F2(-n, -n - 1, ^ + 1; 2, -n - £ - 1; 1) 



(-1) 



^ 4(2n+3)(n+l)! 



(l+2Q!(f+rt/2)! 



[(n/2)!]^ {e+i){e-n/2-iy.{i+2e+ny. 



n even 



,^21-LL 2(n+l)(2n+3)(n+2)! (l+2i')!(l/2+i'+«/2)! , , 

[(n/2+l/2)!p (^+l){^-n/2-l/2)!(2+2^+n)! U'^*^ 

n/2 4(2n+3)(n+l)! 

' -i^ — ...^ ^ n even 



( ■ 2"[(n/2).j 

^ Order (1/^) ^ n odd 



oo 



(5.6) 



(5.7) 



where B{a,b) is the Beta function and we used the Stirling's formula zl ~ \'27rz^^2e ^ in 
the last equation. 



*This form can be obtained by using the decomposition c^,o.o = /n5n,o.o/"^hj which is briefly discussed 
in the appendix. 
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In the case of ni = n2 —> oo, becoming independent of £, is given by^ 

cl,oo,oo = Co'n' 22niff+i)! 3i^2 (-n, -n - 1, | 2, ^ - n; 1) 



^ I (2n+3)[(n+l)!] ^ 
[{n/2) 



+ n:ira-:ir/„ /nMi4 " even 



(«+l)(2n+3)(n+l)!(n+2)! ^ 



(5.8) 



2^"[(n/2+l/2)!] 



^ (-1) — n^, n ^ oo 
vr 

We can use this information to get some measures of the size of the scalar meson. First, 
we can get an exact answer for {£,n) = (oo,0), 

oo „£=oo 
,2/ 2\ _ . "^"AO 



ml{r})e-.oo = ^Yl 



' 22fc(A:!)2(2fc + l) ' 22(fc+i)[(A: + l)!]2 

= sinh'^ (1) - [log(l + \/2) - log 2] = log 2 

Next, (r/)^ 

_>oo can be estimated well enough numerically by considering the first few terms 
in the summation because c^^^ alternating sign and the magnitude of c^^^l'^n de- 
creases for large n: for large even n, 

^ ^ (_l)n/2 2 ^_3/2 _ 

V27r 



This gives 



For general 



^ oo g^=00 

mfe(r/)£_,oo = o y] - 0-697 . 

2 ^ [mnY 



21-1 . 

vr v-^ c. 



ra,0,0 



2 ^ (m„)3 

n=0 ^ 

which increases with i but is bounded by {rf)i^ao- We can also estimate {rf)n^oo by using 
ci,oo,oo for m = n2 = oo: 

mh{rf)n-,oo ~ 0.61 . (5.9) 

More exact results on {r'j)i can be obtained by differentiating Eq. ( p| ) with respect to 
g2 first and then integrating. For the scalar meson with {i,n = 0), 



'Here we used the asymptotic representation of the Jacobi polynomial, 

P^-«(cos^) = cos([n + (a + /?+l)/2]^-W2 + l/4).) ^ ^ 3 ^ ^ ^ ^ ^ 



7rn ( sm ^ ^ 
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where -fff = l + l/2 + l/3 + -- - + l/^is the i-th harmonic number. For example, {r'j)£=i = 
5/12m^ and the squared radius is bounded by 

{rjh^oo = ^ . (5.11) 

The effect of large n is to make (r'j) decrease to some extent, but (r'j) is bounded from 
below : 

, 2(1 - log 2) 

r? ™ = ^ 2^ . 5.12 



Examples: 



In the case of (1, 0| J^(g)|l, 0), 



6 



+ 



6 



12 



6ml 



TT 



2mox 



-12m^ log (m/ix) 



X 



X ^ 



oo 



The logarithmic divergence for x — > stems from the leading behavior for large q^. 
More generally, 



0,712 



m = (-1) 



n2+l 



12 



2n2 +3 



n2(n2 + 2) 2n2 + 3 (n2 + l)(n2 + 3) 



g2 + ml^_-^ g2 _^ ^2^ 



g2 + fhf^^+i 



In the limit of large q'^, the leading behavior of this form factor goes like 1/q'^. 
5.1.2 SO (4) current 

The 50(4) current is associated with the Killing vectors which generate the isometry of S^. 
For simplicity we limit ourselves to the two U{1) subgroups in the maximal torus of S'0(4) 
under which our hadrons are eigenstates, namely the diagonal U{1) subgroups of SU{2)^ 
and SU{2)r. 

In this case, the setup is almost identical to the one in [pO| . We consider the interaction 
of the hadron modes with the canonically normalized 50(4) mode on the D7 brane, given 



in section 2.2, through the vertex 

where iQ is the eigenvalue from acting with the Killing vector v'^da on the mode. 
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When the non-normahzable mode with canonical normalization is put in this vertex 
along with the normalizable modes, we have the matrix element 



{i,n2;p'\J!:\i,nr,p) = Qe^{p + pTiv^^u - q,qu/q'){27T)'6^ |^p^j F^^^i^if) , (5.13) 
F^Z.if) =^ t% qV^K,{qV^)ct>i^^{v)ct>i^^{v) . (5.14) 



Note that ( 5.14 ) approaches the normalization integral for the normalizable modes in the 
— > limit. This guarantees that -Fii^n2(0) = ^n-^,n2■ 
Since (f)\^{v) is a polynomial in ^/v, the exact form of the integral can be obtained 
for arbitrary modes in principle, but the general form is hard to obtain. However, we can 
obtain some interesting results which apply to all or many of the states. First of all, as 



in the flavor case, the form factor goes like 1/(92^^0-1. From Eq. (|1|), q^Ki{q^) 



falls to zero very rapidly and the integration gets the contribution mostly from the region 
< w ^ 1/^2. It is therefore useful to convert the integration variable to q'^v. Now, suppose 
q^ ^ mf^ is large {q^ S> 1). Then the term with the lowest power of q^/v yields the leading 
contribution. The calculation only with the term of the minimal power in v is easy and it 
gives us exactly Eq. (^).^^ 

To compute the two-dimensional Fourier transformation of the form factor is straight- 
forward. 

Again, we can extract the important information without getting into too much computa- 
tional detail. First of all, instead of the previous equation, we consider the following 



2 Jo w + x- 

which has the following relationship with the Fourier transformation we want: 



F(°)(x2) = ^ / ^'tl^^ (5.16) 



F(i)(x2) = 2Ax2f(0)(x2) . (5.17) 
dx^ 



This is proven in the appendix, as are the more general set of relations ( [4.1D and (4.2) 
Now, by partial integration, we have 



i^S.(x^) = ni(x2)iog(i + -i)+n2(x2) 



^"This is not a surprise. In the large limit, the form factor is not affected by the quark masses and is 
governed by conformal invariance. In other words, the equation for the flavor mode can be written as 

d(q^v)^ 4 d{q^v) \ d{q^v) 

Therefore, we can consistently expand G{q,v) — G(,{qy/v) + Gi{q^/v) /q^ + ■ ■ ■ where Go is precisely the 



conformal mode (2.9) 
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y2 



In the small-x region, 112 (x^) dominates because the lowest order term of ni(x^) is ~ 
(x^)^~^, except for the case i = \ for which there is a logarithmic singularity as x ^ 0. 
This behavior is the same as in the flavor form factors. Once again, this is not a coincidence 
since it is related to the behavior of the form factors at large g^, which should be the same 
for both cases by conformal symmetry. 

The large-x behavior is obtained by approximating {v + x^) in Eq. ( |5.15| ). The 

coefficient of the leading 1/x^ term turns out to be non-zero only for n\ = n2 and ni = n2±l 
because of the orthogonality and the recurrence relations of Jacobi polynomials. 



We also compute the 50(4) charge radius squared (rf) of the mesons as given in (3.4). 
In the purely conformal case A = 0, we find a logarithmic divergence. This is an infrared 
effect, stemming from the continuous spectrum of the conformal field theory in the = 4 
sector. The introduction of a nonzero A regulates the divergence. The non-normalizable 
modes for A are now those of Eq. ( p. 161) , and so we obtain 



d 



.-^/^|i^i(.-V^)+/i(.-V^)^#^| 



d^'^'' 'r'''-''''^ ''''''' ""'l^i^^^i 



2m|l^°nx)-2^°^"+4 



Thus (^rfj can be written in closed form for the ground state meson ni = n2 = for any i. 

Note that this has the earlier-noted logarithmic divergence in the limit A — > with rrih 
fixed, but goes to zero in the limit mh — > oo with A fixed. 
For large i, rii = n2 = 0, we find 

which shows the same divergence as before. As rii = n2 becomes large, i fixed, 

(rl) i^|log(!^)+log2 + || , (5.20) 

which also shows the same logarithmic divergence. 

However, the logarithmic divergence is misleading. It is merely an indication that (^rf^ 
is not a good measure in this theory. Since q^/vKl{q^/v) = 1 + g^f log [cp'v)/A + ■ ■ ■ when 
Qy/v <C 1, the two dimensional Fourier transform of the second term is 

/2v 
d'^qe^'^'^q'^vlog {q'^v)/4 ~ — + (UV sensitive terms) . 

This 1/x^ tail leads to a logarithmic divergence, since 

If 11 

(^?) 2 / • — 2 log("l/,/A) . 

J x^ 
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On the other hand, there are other measures of the size which don't suffer from this 
divergence. For example, from the previous equation it is clear that (ri) is finite. Using 



dx ■ 



(x2 + u)2 ' 

we can show that for n = and general i, 

^ ' 2m/, r(<)r(2< + 1) 

which has a limit 

in) TT/{2V2mh) {£^oo) . (5.22) 
This result is not significantly modified for the excited states: 

(ri) — . (5.23) 

ni-+oo mfi 

Examples: 

We begin with the form factor of the matrix element {i, ni = Q\J^^\l, 71,2) in the position 
space. As explained before, the computation of F'^^^ is much easier in many cases and the 
higher spin form factors can be derived from it. We easily obtain 

2' Jo 



Jo \dv J u + x' 



.2 



)n2 f-1 
Jo 



ydx^ J Jq v + x^ 

Here we discover an identity, 

which is derived using a relation proven in the appendix. The rest of the computation is 
str aightfowar d , 

ni = (l + xY+i(-x7-\ %=V+-i^'), 
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As an example, the diagonal form factor for the lowest lying state is 



Fo7=\x2) = -3(5 + 6x2) ^ g " [^2(1 ^ ^2)2] i^g(i ^ 



d 



dx2 



Note that F^^q^ has an interesting limit as i goes to infinity. Since the minimal power 
of the polynomial multiplying the log is at least v^~^, the log part makes the subleading 
contribution in that limit. This is the same for the second summation in ( ^.25 ). The 
coefficient in the first term multiplied by (C/q)^ simplified greatly (^ — 1 — — > 
2^^^^. Hence, the limit is 



^0,0 (X) 



j=Q 



+ x^ 



;i+x2)2- 



(5.26) 



This result can more easily be obtained by noting that the wavefunction-squared in the 
large i limit is peaked at u = 1/2 and becomes 6{v — 1/2). 

For generic i, the diagonal form factor approaches the following form at large ni, 



ni.ni \ J 



1 



ni-^oo x(l + X^y^/'^ ' 



(5.27) 



which can be obtained by using the same approximation used in deriving Eq. (5.12). For 
finite rii there is no actual divergence at x = 0, where the integrals can be done explicitly: 



pWi (0) 

ni,ni v"/ 



2(2ni + 2^ + 1) 
Z^l 



(5.28) 



The computation in momentum space is not as tractable as that in position space. For 
'on2~^(^^)' ^ integration is done in a simple way. We first use an identity 



q^Kjyiqx) 
to write ( ^.14 ) in the following form 



^ ni,n2 



dm 



m? L 

(p + TV? 2 



2 [-1 



dm ■ 



dv 



q2 _|_ ^2 



— y/vJl (m Vw)(/>i_„^ {v)(t>in2 (^) 
^ 



(5.29) 



In the case of £ = 1 and rii = 0, we use partial integration and the recurrence relations of 
the Bessel function in order, to finally obtain 



{1)1=1, 2^ 
^0,712 (9 J 



dm ■ 



m 



+ m 



2 c;„c;„,(-ir('i2+2)(n2 + i) 



m J 



m 



— n2J2n2+2.{m) + J2n2+4:{m) 



(5.30) 



The m integration can be done in the following way. Since 

J,(z) = Mi/«(z)-e'-/7«(-z) 
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we can convert ( |5.3C1| ) to the integration along the contour shown in Fig. and we finally 
have 



n2 + 2 

2 



{-2n2K2n2+3iq) +qK2n2+4{q)} 



(5.31) 



n2+2 



(<!) = E 



(_l)n2-J(„2+j)!(i + l) /2 



j=0 



(ri2-J + 2)! 

In particular, the diagonal form factor for the {i, 111^2) - 

-2^ 12 



2i 



(1,0) 



0,0 



^ 16 192 ^ , , 
l-^ + ^-4i^4(g-) 



(5.32) 




Figure 1: The coutour used in the integration of 

Eq. (Israel) . 



Note that s^4^^~^(g) exactly cancels the 
singular pieces in the expansion of the Bessel 
function around q = 0. This is required in 
order that -Fo,n2 ^o,n2 as q'^ — > 0. Also 
note that the form factor goes like 1/g^ at 
large q, as we have noticed before. 

5.1.3 Energy-momentum tensor 

We can evaluate the matrix element for the 
energy-momentum tensor which couples to 
the boundary metric. The calculation is 
similar to the 50(4) case except that we 
deal with the bulk graviton. We only sum- 
marize our results here. The traceless part 
of the matrix element is given by 



{£,n2;p\T^'''\£,ni;p') = (tj' 



-.2 I P<^ 



per 



-iq' 



ni ,n2 



dv 



{q'^vK2iqVv)} (l)l^^{v)(l)l^^iv) . 



(5.33) 
(5.34) 

(5.35) 



4 Jo 

The analysis is parallel to the SO (4) case. In fact, by recurrence relations between the 
Bessel functions, we have the relation 



(5.36) 



-28- 



as follows from relations proven in the appendix. Hence, the large behavior is dominated 
by the same leading power as the flavor and the 50(4) cases, 

The Fourier transform of the form factor is also easily carried out. We obtain 

2\3 ■ ^ini i^)(t>ln2 i^) , (5-37) 



2 ./ '^'' {v+^^r 



and the Fourier transform also has the corresponding relationship with the 50(4) case, 



1 d 
dx^ 



x^F(i)(x2 



dx2 



X 



^FW(x2 



(5.38) 



also following from results given in the appendix. We can follow the same series of arguments 
to extract the general behavior. It's the sum of a rational function in x and another rational 
one multiplied by log(l + l/x^). The leading power of the log part in the small-x region is 
x2(^~^) again. In the large x region, F^x^) ~ 1/x^, although the coefficient of the 1/x^ tail 
is non-zero only for ni = n2, n2 =b 1, n2 ± 2. 

Because the tail of -F^^) f^Hg ^ow as 1/x^, the quantity (r^) is finite, unlike (rf), even 
as A ^ 0: 

(^ + 1) 



1 

8^ 



For large I we have 



(4) 



1 



1 



6{2e + 3) 




The study of the large ni = ?i2 limit for i = 1 gives 

1 f 1 




1 



-)■ 

mhj 



(5.39) 



(5.40) 



(5.41) 



From the momentum space point of view, the absence of the divergence is related to the 
corresponding absence of the log (f in the small q expansion of the non-normalizible mode. 

Indeed, we have q^vK2{q^/v) /2 = 1 — q^v/'i — q^v"^ log(^2i;)/32 -\ . 

Again we may calculate (r2) for each case that we studied in the 50(4) context. The 
actual computation is unnecessary since we have, from Eq. ( |5.38| ), 



1 /• d^x 

rrih J 2tt 



1 d 
x^ dx^ 



x4f(i)(x) 



1 



2mh. 



dxx2F(i)(x) 



(ri) 



(5.42) 



Hence, we end the discussion by referring to the results of the previous section, in particular 
Eqs. (|]2lD and (|533|) . 



Examples: 



Since F^"^^ can be obtained from F^^^ using Eq. ( ^.38 ), we can use our earlier results on 
the latter to find, for instance, 

2 



p(2K=l 



(x) = -36(1 + 2x^) + 6 



i°g ( 1 + ;^ 
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This is consistent with our prevous analysis on the general cases. In series expansion near 
X ~ oo, we discover that the leading term is 1/x^ and there is a logarithmic divergence at 
X = 0, the latter being absent for i > 1. 

In the momentum space, -Ponf~^(9^) ^ similar expression as the 50(4) case, 

2 \ ^ " 1 

S^nH'l) - - {(4^2(^2 - 1) + q^) K2n^+M + 8qK2n,+M} 



s. 



(2).-^ "^' (-l)"^-^-(n2+j)!(j + 2)(j + l) f2y^ 
"^^^^ fr", (n2-J+2)! {qj 



In particular, the diagonal matrix element of the lowest lying mode yields the form factor 



p{2)^=i. 2n _ 24 
^0,0 ) - 



24 384 

I- -2+ — - qMq) 

q^ 



Note that the term, which is leading at large g^, is consistent with our analysis in 
general. Near q"^ = 0, it can be easily checked that Sn}{q) cancels off the all of the singular 
terms coming from the Bessel function part. Therefore, the form factor is regular at q^ = 0. 



5.2 Spin one hadrons 

We now turn to the spin-one hadrons, type // in Ref. |1C], which are created by acting on 
the vacuum with the operators (Q^^^Q ~ Q^^Q'^)ee- 

5.2.1 Flavor current 

The coupling of the flavor current to two spin-one hadrons descends from the following term 
in the Born-Infeld action. 

58 j d^^gg^^'^g-'^r'^dy^Al^A^Al. (5.43) 

In general, the current matrix element {p\ C'\Jp\p, Cl of a spin 1 meson can be arranged 
into the sum of electric, magnetic, and quadrupole form factors, as in Eq. (^^). Plugging 
the non-normalizable mode for Aa- with no derivative and the normalizable modes for the 
other ^'s in Eq. ( 5.43| ), we get the part of the matrix element. On the other hand, we get 



the Fm part of the matrix element by plugging the non-normalizable mode for Ay on which 
acts and the normalizable modes for the other A's. However, there is no quadruple form 
factor: Fq = 0. Moreover, the electric and magnetic form factors are equal, F^ = Fm = F^f\ 
and they are of the following form: 
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Since Ff>{q^ — > 0) = 1, the static magnetic (/i) and quadrupole {D) moments of the vector 
meson have no anomalous component: 

H = 1 + F„(0) = 2, 

D = ^[F^(0)-2F,(0)] = 



The general form of this form factor is the same as that of the spin-zero case, Eq. (5.2), 
except that the upper limit of the summation is now ni + n2 + 2i + 1 with ^ > 0. Following 
the same argument as that of the spin-zero case, we get the large q'^ behavior of the form 
factor for the general case. In analogy to Eq. (|5.3| ), we find the large behavior is given 

by 

I'^ax ( -,-.^+712'^ fill fill ^ nr ^ (^AA') 

where Ai = £ -|- 3 is the conformal dimension of the spin-1 mode. 

As in 1 24 1, we can classify the form factor into three parts according to the polarizations 
of the incoming and the outgoing hadron states in the Breit frame (the frame in which the 
initial and final hadron have equal and opposite momentum vectors.) In particular, when 
both of the polarizations are longitudinal, the form factor is 

2m'^ ™ ~'~ ( ~'~ Ani"^ ) 

where m is the mass of the hadron. From = and Fq = 0, it follows that Fll ~ 
l/(g2)^^~2. Equivalently, Fll ~ l/{q'^y'^~^ , where ri = Ai — 1 is the lowest twist among all 
operators which can create this spin one hadron. This is consistent with the parton counting 
rule which applies at weak coupling. In this regime, it is expected that Fn ~ 
where p is the number of valence partons; but ri = p at weak coupling. 

Meanwhile Fit and Ftt, where T stands for "transversely-polarized" hadrons, are 
expected at weak coupling to be suppressed each by l/g and l/g^, due to the breaking of 
helicity conservation |2^, |2^ . We similarly, though trivially, find the same behavior at large 
A: Ftt = Fe ^ and Flt ~ qFe ~ 1/ {q^Y^'^''^ . 

It is straightforward to compute the Fourier transformation of the form factor. The 
result is almost the same form as that of the spin-zero case, Eq. (l5.4|) , except that the 
upper limit of the summation is now ni + n2 + 2i + l with ^ > and there is no logarithmic 
divergence for any ^ as x ^ 0. 

The expression for c^^ q q is given by 

cl,o,o = C^!nCli^B{£ + 3,n + e + 2YF2{-n, -n - 1,£ + 3;2,-n - £ - 1; 1) 

4(2«+3)(n+l)! 2(2<;+3)!(2+^+n/2)! „ pvpti 

V ^W2TP (i+r)(e~n/2)U4+2e~\-n)\ At even 



( ^V-^ 2(n+l)(2n+3)(n+2)! (2m)!(3/2+f+n/2)! . , 

V ^) [{n/2+l/2V.y^ (l4-l)(£-n/2+l/2)\(4,~\-2i+nV. ^ 

(-l)n/2 4(2n+3)(n+l)! 

V 2"[(n/2)!]2 ^ 6ven 



(5.45) 



[(n/2)!]2 {l+l){£~n/2y.{4:+2£+n)\ 

fl)(2n+3)(n+2)! (2m)!(3/24 

[(n/2+1/2)!]^ (/+l)(^-n/2+l/2)!{4+2£+n)! 

00 (5.46) 



Order(l/^)^0 n odd 
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In the case of ni = n2 ^ oo, we get the same expression for oo oo the spin-zero case, 
Eq. Q. 

We know that {r'j)e^oo as well as {rf)e^oo is the same as that of spin-zero case because 

''n,0,0 



for £ ^ oo approach the same limit. For any i, 



21+1 f 
n=0 

which has the same behavior with respect to i as in the spin-zero case. {rf)n^oo is also the 
same as in the spin-zero case. 

The mean squared radius of the ground state (n = 0) vector meson of a general i is 

This differs slightly from Eq. ( |5.10| ) but has the same general behavior, and it has the same 
large-^ and large-n limits as Eqs. ( |5.11| ) and ( p.l2| ). 

Examples: 



In the case of (0, 0| J^(g)|0, 0), 



p(/)^=0/-2N _ 1^ _ 192 

^00 W ; - -2 , -2 .2,^2 pj4 



F^P'^^i^) = 12ml [Koimo^) " i^o(mix)] 



^6m|log3, x^O 

Since there is no 1/q'^ behavior for large q'^, FqP^~^ (x^) — > constant as x — > 0. 
More generally. 



'24(n2 + l)(n2 + 2) 



n2 2n2 + 3 _|_ n2 + 3 



V (2n2 + 3) 
5.2.2 SO {4) current 

One part of the matrix element for the SO (4) current descends from the term 



K 

R 



K 

R 



/ • (5.48) 



^^This can be understood without calculating c^~o^ explicitly. Both 3F2(— n, — n — 1, + 1; 2, — n — 
and 3-F2(— n, —n — 1,£ + 3;2, —n — i' — 1; 1) approach the same limit 2^1 [—n, —n — 1; 2; —1) as £ 00. The 
ratio of Beta functions together with the normalization constants also becomes 1 as i ^ 00. 
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Unlike the spin zero case, there is another contribution from 



K 

2R 



R 



Q ^d,A,{{An}{Anr-{A''UAn}} . (5.49) 



The matrix element is of the form given in ( |3.3| ), but the form of ( ^.48| ) and ( ^.49 ) implies 
once again Fg = = F^^^ and Fq = 0. The general form is similar to Eq. (|5.14| ), 



^ Jo V 

Following the same argument as the scalar meson case, we rediscover Eq. ( ^.441) . 
As before, the Fourier transform is obtained by ( 5.17 ) and 



(5.50) 



J.II All 



(5.51) 



ni(x^)iog 1 + ^ +n2(x^) 



x^ 



ni(x^) 



It is easily seen that the general properties of the vector mesons are the same as the scalar 
mesons, though the degree of the polynomials is different. Less obvious is the fact that 
both at large ^, for ni = 71-2 = 0, and at large ni = n2, for fixed the spin-one form factors 
approach the same limits as the spin-zero form factors, Eqs. ( |5.26| ) and ( ^.27 ).-'^^ 

The charge radius {^rf^ is similar in form to that found for the spin-zero case; it also 
has an infrared divergence in the A — > limit. 



For large d. we find 



{rl) 



mt 



V A ; 



+ ■7,{H2i+A - He+2) + 7 



£^00 m 



(5.52) 



(5.53) 



Note the large i limit is the same as the spin zero case, Eq. ( 5.19| ). This follows from the 
above-mentioned identity of the form factors in this regime. Similarly, we find that the 
large ni limit for generic fixed i is identical to that of the spin zero case, Eq. ( p. 20 ). 
Meanwhile, for the ni = n2 = states, we have 



in) 



TT r(£ + |)r(2^ + 4) 



(5.54) 



^Eq. (|5.2^) is however replaced with F^^Wi^ (0) = 2(2ni + 2£ + 3))/(^ + 1) 



-33- 



Again, the identity of the spin-zero and spin-one form factors at large i implies (ri) 
{7r/2\/2)mJ^^ , as in the spin zero case Eq. (5.22). For fixed i, the large-ni limit is l/nih, as 



for the spin zero result Eq. (5.23) 



Examples: 



Following the same line of computation as in the spin zero case, we find the identity 



and a similar expression for the case of the ground states 



(0)£+n2/ 2 



piu^i;/ 2\ _ fr'II\2 
-^0,0 (X ) - (C^o ) 



ni(x^)iog 1 + ^ +n2(x^) 



X^ 



ni(x2) = (1 + x2)^+i(-x2)^+i, n2(x2) = P++(x2) 

where P++ is defined in Eq. (l5.25| ). In particular, 



5(1 - 8x2 _ gg^4 _ gQ^6^ _^ [x^(l + x2)2] log ^1 + -^2^ 



Note that this is the lowest state that the form factor calculation is meaningful. Since the 
spin one hadron is in the {i/2,i/2) representation, the £ = hadrons are neutral under 
5*0(4) transformation. 

In the momentum space, we compute the matrix element {£ = 1, n2| J4 |1, 0). The form 
factor is 



<r(^") = 3./5(2n2 + 5)^"^ + ^ 



(^) - 3i?2n2+7(4/g') ( ^ ) K2n2+7iq) + R2n2+8i^/ f)K2n2+8iq) 



R2n2+7iz) = (n2 + 2)^2 + u' ^2 + 3 



3 . + 16(2n2 + 7)r^ + ' 



Here Sn2{Q) is minus the singular part of the rest of the expression. In particular when 
712 = 0, 



(5.55) 



p(i)£=i._2. _ 23040 

^00 il ) - -^6~ 



96 3840 21g + g3 



24g2 + g4 
48 



We have a leading term at large q'^, as predicted. 
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5.2.3 Energy-momentum tensor for the spin one case 



We can also compute the energy- momentum tensor matrix element for the spin one hadrons. 
It comes with the similar tensor structure as ( |5.33| ) and ( ^.341 ), and the form factor can be 
written in the same way as ( ^.351 ) except that the metric factor is now different. As before, 
we can derive these form factors from the 5'0(4) form factors using Eqs. (|]3|) and ( ^l38|) . 

The general properties are similar to the spin zero case. The form factor has the same 
large behavior as the SO(4) current, l/g^^^"*"^). For n = and general 



{rl) 



8ml 



8mi 



(^ + 3) (AV 
6{2£ + 5) \mhj 



^ 12 



rrih 



There is no logarithmic divergence. Again, since the spin-one and spin-zero form factors 
are the same in the £ — > oo limit (ni = n2 = 0) and in the large rii = n2 limit for fixed £, 
our results for (r^) agree with the spin- zero case in these computations, namely Eqs. ( ^.40 ) 
and ( 5.41| ). Moreover, it is again true from Eq. ( ^.42 ) that (r2) = (ri)/2, and so for ni = n2 
and general i, (r2) can be obtained from Eq. ( |5.54 ). The large i and large rii = n2 limits 
can be similarly read off from Eqs. ( ^.22 ) and ( 5.23 ). 



Examples: 



Using Eqs. ( ^.36| ) and ( 5.38 ), we find 



Fr 



(2)£=0._2n 



0,0 

92160 



576 

;74 



32 n-^ rr 



120 5760 48g2 + g4 48f + f . 



and their Fourier transforms are 



^ir°(x^) = + 12x^ + 12(x)^) - 12 (^)' [x^(l +x^)] log (l + i,) 



F^t\^') = 10(1 - 12x2 _ 150^4 _ 180^6) + go [^8(1 + ^2)2] ^ J_ 



6. Implications 



Combining the results of |10| with those of the previous sections reveals a number of un- 
familiar patterns. Taken together, they confirm that this is a class of bound states quite 
unlike any previously studied. 



6.1 The quarkonium spectrum 



Of course the first surprise involves the mass spectrum itself |1C]. As we noted earlier, we 
would not normally expect quarkonium of spin < 1 to be so much lighter than quarkonium 
states with higher spin, and certainly not to have a mass of order ruh ~ A~^/^mQ. Also 
surprising is that the masses of states with radial quantum number n, and with I of the 
light $i particles added to the (s)quark and anti(s)quark, is approximately linear in {n + £). 
One might wonder if a constituent quark model might apply to these hadrons, though this 
would require quark constituent masses of order ruh (much less than the bare quark masses) 
and constituent masses for the of order mh (much greater than any confinement scale A, 
if there is any low-energy confinement at all.) 

All of these facts are of course straightforward to understand from supergravity. In 
the supergravity limit, the scale nih is the only one which appears in the equations solved 
by |10|. Moreover, the n + i dependence is at least partially guaranteed by the extra- 
dimensional Kaluza-Klein-like structure. The quantum number n sets the number of nodes 
of the hadron's wave function in the r direction, while £ sets the number of nodes on the S^, 
so linear dependence in n when I and in £ when £ 3> n is natural. But these arguments 
give no insight into how to derive these facts independently from quantum field theory. 

Similarly, the large binding energy of these states is no surprise from string theory. The 
(s)quarks are strings connecting the D3 branes to the D7 brane, and in the supergravity 
limit are strings extending from the D7-brane to the horizon of AdS. They are much longer, 
and indeed more massive, than the 7-7 strings, which are massless on the D7 brane. But 
from the field theory point of view, this is still completely mysterious. The strange nature 
of this binding is highlighted by the fact that it is vastly reduced for quarkonium involving 
two (s)quarks of significantly different mass, or, even more remarkably, between (s)quarks 
whose mass parameters mi and m2 differ only by a phase. The masses of such states are of 
order \mi — m2|. The phases of the mass parameters appear in the interactions of the field 
theory, but how they conspire to make deep binding possible when mi = m2 is unknown. 



6.2 The size of quarkonium 

Perhaps our most striking new observation is this: the form factors we computed indicate 
that these states all have sizes of order = y/X/mq, for all n and £. 

From the supergravity point of view, this is not that hard to understand. Just as the 
supergravity equations ensure that nih is the only scale which can determine the masses of 
hadrons, so it is the only scale which can appear in the form factors in the nonconfining 
theories (where A = 0.) Thus for A = it must be that F{q^) is really F{(f') (recall 
q = q/rrih) and F{x±) is really -F(x) (recall x = mhX±.) We therefore should not be 
surprised that mq does not appear in any of our expressions, since any such appearance 
would require that a' appear in the bulk physics, through some string theoretic effect beyond 
supergravity. 

But nih is not the only possible scale. A hadron's size, especially if it contains at least 
one $ particle, certainly could be infinite in the A — > limit. We might have expected that 
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at least some measure of the charge radius of these hadrons would have been of order A~ . 
This does not happen (at least for the conserved currents.) 

What is strange about this result is that it differs greatly from what we would have 
expected at small 't Hooft coupling. Let us consider pure (s)quarkonium first, and then 
(s)quarkonium with $ particles added; at small 't Hooft coupling these are very different 
systems. 

6.2.1 Pure (s)quarkonium 

The pure (s)quarkonium system, in which the (s)quark and anti(s)quark are in a color-singlet 
state, is hydrogenic — or more precisely, positronium-like — at small 't Hooft coupling. At 
A = there is an exactly Coulombic potential between the (s)quarks, so at small A the size 
of (s)quarkonium is certainly ~ (XniQ)"^. As A increases, the system becomes relativistic 
(for small angular momentum), so this estimate breaks down. We do know on general 
grounds that the size of low- lying states will be r ~ fW/mq, where / is an unknown 
function that behaves as 1/A at small A. At large A, /(A) could certainly be of order 1, or 
even smaller; we have no preconceived notion from the field theory point of view of how it 
should behave. 

In principle it might have been the case that the quarkonium bound state approached in 
appearance a point-particle, with size much less than its inverse mass, in the large A limit. If 
/ continued to shrink, or even went to a constant, this would have been the case. However, 
this does not happen; /(A) reaches a minimum somewhere around A ~ 1 and then begins 
to grow. The (s)quarkonium states remain as large as, or larger than (for highly excited 
states), their inverse masses. In this sense they retain the fluffy properties of composite 
objects even at large A. 

6.2.2 Generalized (s) quarkonium 

The behavior of the states containing $ particles is even more difficult to understand. 
For example, consider the case with one $ added to Q and Q. In this case the (s)quark 
and anti(s)quark arc combined in the adjoint representation of color. We might naively 
expect, therefore, that they repel, as they do at leading order in A; however the repulsion 
is suppressed at large N, so instead we should think of them as noninteracting. The only 
interaction between them is induced by the light $ particle. In this sense, this system is 
like a hydrogen molecule, but without the repulsion between the protons. 

At small 't Hooft coupling it is straightforward to carry out a Born-Oppenheimer com- 
putation of this object. As is easily seen, the light # particle has a wave function that 
spreads out over a distance scale L$ ~ 1/Am$. This answer is nearly independent of the 
distance between the Q and Q; it is manifestly true both when the Q and Q are well- 
separated from one another (in which case the $ can be in one of two patches of size -Zj$, 
one near Q and one near Q) and when the Q and Q are placed at the same point in space 
(in which case the $ has a hydrogenic wave function of size ^ L^.j 

The two heavy particles now move in an effective potential induced by the fact that 
both are attracted to the ^ particle. It is easy to show that the average size Lq of the Q and 
Q wave functions is large compared to l/Amg and small compared to 1/Am#, independent 
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of the precise details of the potential. In the limit m$ — > 0, uiq and A fixed, both L$ and 
Lq diverge. This also happens in the limit A — > with fixed Q and $ masses. On the other 
hand, in the limit mq — > oo, m$ and A fixed, L<j, is constant and Lq — > 0. These results 
would be refiected in the size of the hadron as measured by the fiavor current (which would 
see only the heavy particles) and by the S'0(4) current (which would be sensitive to <I> as 
weh.) 

What we have learned about the large-A regime is very surprising. We have found 
that both the fiavor and SO (A) currents see a hadron whose size is of order ^/X/uiq. This 
is parametrically smaller than either current would see in the small A regime. Taking 
interesting limits makes this especially clear. For m$ 0, mq and A fixed, the size of the 
hadron is fixed from the point of view of both currents. For mq oo, m$ and A fixed, both 
currents see a hadron shriking to a point; and this is even true when mg — > oo, — > 0, 
with, say, rnqm^ fixed. This is a striking phenomenon not seen previously in quantum 
field theory, to our knowledge. The light, or even massless, $ is somehow trapped by its 
interactions with the heavy particles at a size scale of order mq. 

Comparing this behavior with the small-A regime, we see that from the point of view of 
SO(4) (and, to a lesser extent, fiavor) the derivative of the hadron's size with respect to A 
is extremely large, and negative, near A ~ 1. As m$ and/or mq oo, the dependence 
of the size and shape of the hadron on A is apparently nonanalytic. 

6.3 The meaning of the divergence of {rf) 

Now let us turn to the issue of the logarithmic divergence in (rf). It is conventional wisdom 
that a good measure of size, given a form factor F{q^), is (r^) oc d[F[q^)\/ d{q^) at = 0. 
However, this idea relies crucially on the exponential falloff of all charge distributions in 
position space. In this model, the S'0(4)-charge distributions in position space fall off only 
as a power of radius, because the (s)quarks are coupled to a conformally-invariant sector 
that carries SO (4) charge. No matter how small the coefficient of might be in a 

transverse charge distribution for the SO{4) current, it always leads to a divergent (r^). 
This proves not that the hadron is infinite in size but that (r^) is a bad measure to use. 
Indeed other form factors, even in the = 4 sector, have finite (r^), as does (r) for (5*0(4). 

More physically, if one simply compares graphs of fiavor form factors and of 5*0(4) form 
factors in position space, as in Figs. ^, |^ and|^, one finds they can be difficult to distinguish, 
and are never very different. Instead one sees that all the hadrons have a core of order mj^^ 
in size, and that their power-law tails are only a small part of the structure of the hadron. 
In particular, for large ^ and n = the fraction of SO (4) charge stored inside |x| < 2 is 
8/9; for large n and fixed I the fraction is almost the same, Ijyfh. 

This fact has implications for string theory. It has long been said that strings are 
infinite in size, because of the divergence found in d\F((^)\l d{c^^ at g'^ = for the energy- 
momentum form factor. However, this conclusion might be erroneous. Strings, which couple 
to massless gravitons, ete., might well have charge-distributions and energy-momentum 
distributions with power-law tails. It would be worthwhile to revisit the question of the size 
of strings in the light of this observation. 
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Figure 2: The diagonal flavor and 50(4) form factors, in position space, for the £ ^ 1, n ^ 
spin-zero hadron. 
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Figure 3: The diagonal flavor and S'0(4) form factors, in position space, for the ^ = 10, n = 
spin-zero hadron. 

6.4 The structure of quarkonium 

As we have seen, all of the hadrons in question appear to have a core and a power-law tail. 
The detailed shape of the core, though not its radial width, depends on n and somewhat 
on i. The tail has a fixed power and a normalization which is rather insensitive to n and i. 

The flavor form factor tells us where the (s)quarks are located, and is insensitive to 
the adjoint matter. Its properties indicate that the (s)quarks are located in the core. The 
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Figure 4: The diagonal flavor and 5*0(4) form factors, in position space, for the ^ = 1, n = 10 
spin-zero hadron. 

50(4) current, which is sensitive to the $ fields and gluinos as well as the (s)quarks, sees a 
tail for any L That this is true even for 1 = 1 (spin-zero), for which there are no valence $ 
particles in the hadron, suggests that the tail is populated by a sea of gluons, gluinos and 
<I> pairs. (Results from [^] suggest that the sea dominates hadrons at large A, in contrast 
to expectations at small A, where valence particles play a more significant role.) The lack 
of any dramatic difference between the tails of £ = 1 and (. = 2 states suggests that the 
valence <I> particles do not strongly change the structure of the hadron. 

At large ^ the 50(4) form factor is independent of ^, while the charge Q of the hadron 
is proportional to I. Raising I by one corresponds to adding a <I> particle into the state; that 
this has no effect on the shape of the hadron suggests that for large £ the valence $ particles 
are all coincident in position and distribution, affecting the 50(4) charge distribution only 
through an overall scaling with 

Increasing n does change the shape of the hadron; a somewhat larger fraction of the 
flavor charge, 5*0(4) charge and energy are located in central part of the core. However, 
even in this case the fractions of the charges located in the tail differs little between large 
and small n. It does not seem, then, that large n has much impact on the basic structure 
of the state. 

In both the large-^ and large-n limits, in which the hadron mass also becomes large, 
the form factors have a simple structure. The core size is of order m'^^ (which is much 
larger than the hadron's inverse Compton length), and the tail has a fixed normalization 
and power law. This suggests that the hadron acts like a bag; the ^ particles, trapped 
inside, have a symmetric wave function with radius of order . The puzzle, as before, 
is that the size of this bag is set by 1/mh, rather than the mass-scale of the $ particles 
themselves. 
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The small size of the bag is directly tied to the deep binding of the state, in the following 
sense. Suppose a <^ particle's wave function spreads out to a width d, where 3> d 3> 
rrih- This arrangement suffers a huge energy cost. In the stringy dual description of this 
phenomenon, this corresponds roughly to a part of the 7-7 string which represents the 
hadron descending from the D7 brane, which lies at the AdS radius L = mqa', down to 
an AdS radius of order r = ^/Xa' /d, and then returning to the vicinity of the D7 brane. 
This would have an energy cost not of order 1/mh, as one might expect, but of order 
2{mQ — y/X/d). Thus the same physics which makes the 7-7 string nearly massless also 
provides a mechanism by which it can energetically trap the $ particles. Similarly, if the 
<I> wavefunction remains small but moves a distance d from the center of the bound state, 
this costs of order 2^/\m'^d = 2{mQ/ X)d. This linear growth of the energy with d is a 
less drastic effect than the previous one, but is still substantial, and is consistent with the 



Regge-behavior found in |10| for states of moderate spin, which are of linear size larger 
than l/ruh- As emphasized earlier, we still have no idea from the field theory point of view 
why the deep binding is present, but it is clear that the small size and deep binding are 
correlated. 

6.5 The case of the missing form factors 

Finally we turn our attention to another key observation, which was first made by Son and 



Stephanov |27] in the context of a recent model [1^. These authors noted that there are no 
anomalous magnetic form factors or quadrupole form factors for the spin-one particles in 
their recent model This is in contrast to hadrons in QCD, such as the p meson ||2^. The 
reason for this can be traced to the supergravity limit. The couplings which are required 
are present in the string Lagrangian but are subleading in the a' expansion. The hadrons 
in our model also share this feature. For instance, anomalous flavor magnetic form factors 
Fm — Fe could only be generated by a term a'iiF^F^ F^ in the D7-brane Born-Infeld action. 
This term, if present, could only generate a contribution to F^ — F^ which is suppressed 
by 1/\/A. In a supersymmetric theory, the coefficient of this term vanishes identically. But 
even were supersymmetry broken at order one, so that the F^ term had a coefficient of 
order one, the existence of a low-energy gravity limit would still ensure that the anomalous 
magnetic form factor was extremely small. The quadrupole form factors are suppressed for 
the same reasons. These facts should have general implications for the structure of hadrons 
at large 't Hooft coupling, but we have yet to understand precisely what those might be. 
In particular, despite the fact that these hadrons have a nontrivial size, there is a sense in 
which they retain some of the properties of fundamental point particles in the large A limit. 

6.6 Unansw^ered questions 

There are many additional computations to consider. A key issue in these theories is to 
understand which aspects of hadronic physics are determined by scale invariance in the 
ultraviolet. Such aspects may apply approximately in QCD. We have obtained some results 
which probably follow from conformal invariance, possibly in combination with the large-A 
limit; it remains to prove this connection and understand which conditions are necessary 
and which sufficient. 
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The shapes and sizes of hadrons in this particular model may or may not reflect prop- 
erties of typical light hadrons found in the large A limit; this remains to be explored. Our 
results presumably can be used to help interpret the structure of light hadrons purely in 
the = 4 sector. Previously, the lack of a second scale in theories such as M = 1* has 
limited our ability both to calculate and to interpret form factors. In the present paper, 
we were able to separate the hadron mass scale nih and the conflnement scale A, to great 
advantage. 

On the other hand, the internal structure of these (generalized) quarkonium states still 
is unknown. The methods of ||20[ | can and should be applied to them, in order to better 
understand what role the (s)quark and anti(s)quark are playing. It would be even more 
effective to study the full off- forward distribution amplitudes [16, 17, 18, which can 
combine these pieces of information into one package. 

Finally, there are a large number of still more challenging computations to do that 
go beyond the supergravity approximation. What happens to these (s)quark-anti(s)quark 
bound states when one quark mass is much larger than the other? or even when their 
masses have the same magnitude but different phases, which changes the interactions and 
the binding energies drastically? What is the structure and the dynamics governing the 
higher spin states? Do any of these stringier states in any way resemble their counterparts 
in QCD? Can we understand, perhaps even quantitatively, how states with various spins 
metamorphose from their small-A to large-A forms? As of now, the answers to these ques- 
tions can only be guessed at, but many of them seem tractable, and there is hope that 
considerable progress is possible even with present-day techniques. 
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A. Decomposition of the Form Factors 

A matrix element {b,p'\J^{q)\a,p) can be written in terms of a sum over many hadronic 
states in the supergravity limit. In particular, the following decomposition of the form 
factor can be used: 



as illustrated in Fig. ^. This decomposition of the form factor into a sum over resonances 
is not justified a priori in a general confining field theory. However, in the large A'^, large A 
limit, it is exact, with corrections of order 1/A^ and 1/\/A. 

For the spin-one current matrix elements, the fn are the hadron decay constants of the 
spin-one hadrons created by the current 

{0\J^iq)\n,p,e) = /„e, U'' - ^) S\p - q) , 
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Figure 5: Decomposition of the form factor into a sum over hadron states, as in Eq. ( A.l). 



where |0) is the vacuum of the thoery, and |n,p, e) is the spin-one hadron state with mass 
m„, momentum p and polarization created by the current operator J'^. The fn are 
associated with the Wronksian between the normahzable and nonnormahzable mode.^^ For 
instance, consider the modes created by the flavor current, which satisfy the differential 
equation for (t>^\p) with ^ = 0: 



+ 



m 



(1 







This is already written in a self-adjoint form, so the Wronskian of two solutions is 



W[(j)i,(l)2] = (0102 - 0201; 



p3 



(A.2) 



where / is a constant independent of p. The fact that f = fn when m = rUn follows 
from the form of the two-point function of the current. We may evaluate /„, by evaluating 
the Wronskian at any value of p that is convenient. We do so at /? — > oo, using the 
fact that the nonnormahzable mode (j)non associated to a conserved current goes to a q^- 
independent constant as p — > oc, while tli6 corresponding normalizctble mode 077, goes a^s 
0n(p) = inP^'^ + order(p"3). 



/„ = lim p-^W 

p— >oo 



where 0non 



,„(m„), 0„] = (pnoniP OO) 

(pnonip oo). For the flavor-current case 



P 



dp 



p— »oo 



/„ = {-ir2{LygsR')C^i. 



The gnab are the coupling constants between this hadron and the incoming and outgoing 
scalar hadrons, which are given in almost the same way as the matrix element, except that 



^^We thank D. Son for sharing his derivation of this fact. 
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instead of a nonnormalizable mode of the gauge boson or graviton corresponding to the 
current, with arbitrary momentum q'^ , we plug into the trihnear vertex a normalizable 
mode corresponding to the hadron n, on shell. For example, in the decomposition of the 



form factor Eq. (5A), the corresponding gn,ni,n2 given by 



9n,n\ 



712 



58- 



1 I 



Combining this result and the previous one with Eqs. (|5.2|) and ([A.l|), we can now compute 



'-n, 711 ,712 fn9n,ni,n2/^h' 



B. Some Identities 



We begin with proving Eqs. ( |5T7|) and ( pSj) . It's enough to prove the following identity. 
Let = V + a where a = x^. Then for an integer S > 0, 



1 



(B.l) 



where we have used the fact that v = Q — a and S derivatives of C"^"") 1 < n < S, all vanish. 

Let's take the two-dimensional Fourier transform of this with respect to d^x, converting 
it to a function of (f, where = q^jmh is conjugate to . The Fourier transform of 1/C 
is 



Now we are ready to show that Eq. ( [4.2| ) implies (U), and vice versa. 

^ (-1)^ e 5^ „„ (\\ 1 
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This relation also follows from 
2tt J d{ 



as 1 !• aS 

-2px2^F(x2) = —J d|x|2d0cos(x/^cos0)^^^x2^F(x2) 

1 1 d|x|Vo(vW)^x^^F(x^ 
d\x\'^de 



(-1)^ 



^5 



d\x\'^de 



:2\S 



x^ 



QS 



/ -2\S TPf -2\ 



d{q 



■2\S 



x2'5f(x^ 



where we have used the fact that the integral is even under cos0 — cos 6, and assumed 
the function -F(x^) is sufficiently convergent at x — > and x — > oo to allow the integration 
by parts (which is the case in our applications of this identity.) 
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Another useful relation that we used in derivation of Eq. ( |5.24| ) is 



9" v'^ _ 1 9-^ 5 5 1 _ 1 9-^ 5. a" 1 
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